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Abstract 


Precipitation is the main component of the hydrological cycle, which is very important in many 
studies, including those focusing on management of irrigation and drainage systems, crop 
performance, environmental issues, floods, droughts, etc. Therefore, this study eveluated seven 
GCMs models of Coupled Model Intercomparison Project phase 6 (CMIP6) regarding the 
simulation and projection of precipitation changes in the northwest of Iran under three scenarios of 
SSP in three periods (2021-2050, 2051-2080, 2081-2100) and compared them to the historical 
period (2014-1985). The trend of precipitation changes was calculated using Mann-Kendall test 
and Sen’s slope estimator. The Linear Scaling Bias Correction (LSBC) was used to downscaling 
the GCMs data and the validity of different models was evaluated using RMSE, MAE and R2 
indices. The results of performance evaluation of CMIP6 models showed that among the studied 
models, MPI-ESM1-2-LR model with an average R2 of 0.86 and RMSE equal to 19.7 at the 
regional level were more accurate than other models in the simulation of precipitation. The results 
of projection also showed that precipitation according to SSP1-2.6 scenario will increase by 2.6% 
in the next three periods and on average in the end of 21st century and at most stations it is 
significant at the level of 0.01%. But according to SSP3-7.0 and SSP5-8.5 scenarios, the amount of 
precipitation will decrease by 14.5% and 3.6% in the end of 21st century, respectively, with the 
highest decrease being related to the rainy areas in the southwest of the region. In general, 
according to different scenarios, in most of the studied areas, the precipitation trend is decreasing 
until the end of 21st century, and it is necessary for officials and planners to adopt the required 
strategies to adapt to the resulting climate change. 
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حکیده 


بارش اصلی‌ترین جزء چرخه هیدرولوژی است که در بسیاری از مطالعات ازجمله 
Ca pte‏ سیستم‌های آبیاری و زهکشی. عملکرد محصولات کشاورزی. مطالعات 
زیست محیطی. سیلاب. خشکسالی و غيره از اهمیت زیادی برخوردار است. لذا در اين 
پژوهش به ارزیابی Y‏ مدل GCMs‏ از CMIP6‏ در شبیه‌سازی بارش و پیش‌نگری 
تغییرات بارش در شمال غرب ایران تحت سه سناریوی واداشت اجتماعی-اقتصادی 
(SSP)‏ در سه دوره (۲۰۲۱-۲۰۵۰ ۲۰۵۱-۲۰۸۰ ۲۰۸۱-۲۱۰۰) نسبت به دوره پایه 
(۱۹۸۵-۲۰۱۶) پرداخته شد. روند تغییرات بارش با استفاده از آزمون من کندال و 
تخمینگر شیب سن محاسبه گردید. جهت مقیاس کاهی داده‌های 5 نیز از رویکرد 
مقیاس بندی خطى (LSBC)‏ استفاده شد و صحت سنجی مدل‌های مختلف با استفاده 
از سنجه‌های متداول آماری» MAE RMSE‏ و 82 ارزیابی گردید. نتایج حاصل از 
ارزیابی عملکرد مدل‌های CMIP6‏ نشان داد که از ميان مدل‌های موردبررسی» مدل 
MPI-ESMI2-LR‏ با متوسط R2‏ برابر با ۰/۸۲ و RMSE‏ برابر با ۱۹/۷ در سطح 
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نتایج حاصل از پیش‌نگری تغییرات بارش نيز نشان داد بارش بر اساس سناریوی 
SSP1-2.6‏ در هر سه دوره odul‏ و به‌طور متوسط تا پایان قرن ۲۱ به میزان 7/1 درصد 
افزایش بيدا خواهد کرد و در بیشتر ایستگاه‌ها در سطح ۱ درصد معنی‌دار شده 
است؛ اما بر اساس سناریوهای SSP3-7.0‏ و SSP5-8.5‏ میزان بارش تا پایان قرن ۲۱ به 
ترتیب ۱8/۵ و ۳/۰ درصد کاهش lay‏ خواهد کرد که بیشترین میزان کاهش أن نيز 
مربوط به مناطق پر بارش در جنوب غرب منطقه است. میزان کاهش در مناطق غربی و 
جنوت غربی در سطح v/v)‏ درصد معنی‌دار شده As]‏ درمجموع با توجه به 
است و لازم است مسئولان و برنامه ریزان راهکارهای لازم جهت سازگاری با تغییرات 
اقلیمی حاصل را اتخاذ نمایند. 


کلیدواژه‌ها: بارش, تغيير اقليم» شمال غرب. مقیاس بندی خطی. CMIPOSSPS‏ 


١‏ - مقدمه 


كرمايش جهانى در حال حاضر يذيرفته شده است (ژانگ و همکاران » ۲۰۱۷). یکی از نمودهای كرمايش 
جهانى و تغيير اقليم افزايش فراوانى و شدت دورههاى همراه با دماهاى بالا و يايين است و برخى مطالعات 
نشان دادهاند كه تأثيرات فرينهاى اقليمى آشکارتر و مستقيمتر از ميانكينهاى آب‌وهوایی است (هيأت بين- 
الدول تغيير اقلیم » ۲۰۱۳). دمای بالاتر نگهداشت رطوبت جو را افزايش داده و اين امر به نوبه خود می تواند 
stabil ay ona‏ بارشو فاق کین قوف OL)‏ و OHA Ob aa‏ فك انيف DIR‏ کات و ting‏ 
فضايى رويدادهاى فرين آب‌وهوایی افزايش LL‏ (هورتون و همکاران» ۲۰۱۵). مسئله تغيير اقليم امروزه 
به‌عنوان یکی از شايعترين مباحث علمى و حتى سياسى-اجتماعى مطرح است. اصولاً تغيير و حركت جزء 
وجودى سامانه‌های طبيعى است. انسان در تلاش است كه با ابزارها و دانستنىهاى خود قانونمندی اين 
تغييرات را کشف كند و بدينوسيله اقدام به پیش‌بینی رفتارهای بعدى سامانه آبوهوا نمايد (عزیری,۱۳۸۳). 
فرايند تغییر اقليم بهويزه تغييرات دما و بارش مهمترين بحث مطرح در قلمرو علوم محيطى می‌باشد. اين يديده 
به [Jo‏ ابعاد علمى و كاربردى (اثرات محيطى و اقتصادى- اجتماعى) آن از اهميت فزاینده‌ای برخوردار است» 
جراكه سيستمهاى انسانى وابسته به عناصر اقليمى مانند کشاورزی. صنايع و امثال آن بر مبناى ثبات و پایداری 
اقلیم طراحى شده و عمل می‌نمایند (ذهبيون و همكاران. (YA‏ در اين بين نيز بارش مهمترين متغير 
Zhang et al‏ 1 

2 IPCC 


3 Yang et al 
4 Horton et al 


سال سيزدهم پیش نكرى بارش با استفاده از مدلهاى 081186 تا يايان قرن 7١‏ در dius‏ غرب ايران Wo‏ 


آب‌شناختی است که ييوند ميان جو و فرآیندهای سطحی را برقرار می‌سازد و ناهنجاری‌های مثبت و منفى آن 
سبب بروز سيل و خشكسالى می‌شود (مسعوديان و همكاران. ۱۳۹۳). برای ارزيابى تغييرات اقليمى 
مدل‌های جهانى اقليم (GCMs)‏ طراحى شدهاند. اين مدل‌ها خصوصيات فیزیکی. چرخش‌ها و حركات 
اتمسفرى را تحليل مىنمايند و به دنبال آن متغيرهاى هواشناسی را در شبکه‌های خاص شبيهسازى می HES‏ 
همجنين اين مدلها در حال حاضر قوىترين ابزار به‌منظور توليد سناريوهاى اقليمى می‌باشند (شائمى و 
حبیبی نوخندان. (AYM‏ مدل‌های GCM‏ متعددى توسعه يافته است و به دليل تفاوت ساختار اين مدلها و 
تفاوت شرايط اولیه‌ای که تحت آن اجرا می‌شوند. نتايج آنها حتى برای یک سناريوى انتشار OLS‏ متفاوت 
است. لذا برای استفاده در مطالعات منطقه‌ای» لازم است خروجی اين مدلها با استفاده از یک روش مناسب 
ريز مقياس شود (كى و همكاران » 9 با توجه به اهميت بارش لازم است اين يارامتر در مناطق مختلف 
برآورد شود. پیش‌بینی بارش برای مناطقى كه اقتصاد آن‌ها مبتنى بر كشاورزى است» بسيار مهم است. از 
طرفى با بيش بينى بهموقع بارندكى امكان مقابله با سیلاب» خشک‌سالی و كاهش خسارت ناشی از آن‌ها نيز 
مهيا می‌شود. در اين راستا فريرا و همكاران' (۲۰۱۸) به پیش‌نگری تغييرات بار شهاى تابستانى در جنوب 
شرقى ابالات‌متحده پرداختند. نتايج حاكى از افزايش قابل‌ملاحظه بارش به دليل افزايش دما و شار بخارآب 
وف SS‏ فیکار (P989‏ هیرشن الراك کیرات Clays DU de Los‏ پو تا در عند با Jald‏ 
سناریوهای RCP‏ پرداختند. نتایج OUS‏ داد دما و بارش در دوره آینده افزایش بيدا خواهد کرد. جیانگ و 
همکاران* (۲۰۲۰) به پیش‌نگری تغییرات بارش در آسیای مركزى بر اساس ۱۵ مدل GCM‏ از “CMIP6‏ تحت 
٤‏ سناریوی ۹5۳" پرداختند. نتايج نشان داد بر اساس تمام سناریوها میانگین بارندگی سالانه افزايش خواهد 
cab‏ ژانگ و همکاران (۲۰۲۱) قابلیت‌های ۶۲ مدل از CMIPG‏ را جهت بازتولید دما و بارش در سين 
کیانگ جين را مورد ارزیابی قرار دادند. نتایج OUS‏ داد اين مدل‌ها توانایی خوبی در شبیه‌سازی دما و بارش 
دارند و استفاده از cus‏ جند مدل نيز عملكرد بهترى نسبت به مدلهاى منفرد دارند. مجدى و همکاران" 
(۲۰۲۲) به پیش‌نگری تغييرات دما و بارش در خاورميانه و شمال آفریقا با استفاده از متوسط ۲۳ مدل GCM‏ 
و دو سناريوى SSP‏ يرداختند. نتایج SE‏ از افزايش دما و كاهش بارش در بيشتر مناطق موردمطالعه است. 


صلاحی و همکاران (۱۳۹۵) به پیش‌بینی تغييرات دما و بارش با استفاده از داده‌های خروجى مدل 
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3 تحت سه سناریوی انتشار در حوضه آبریز درياجه پرداختند. نتایج حاكى از کاهش ۸/۲ درصدی 
بارش در دهه ۲۰۵۰ در حوضه موردمطالعه است. خورشیددوست و همکاران (۱۳۹۷) به پیش‌بینی بارش هاى 
سنگین در سواحل جنوبی دریای خزر در دوره ۱۳۸۹ تا ۱۶۰۸ پرداختند. OLS quU‏ داد در دوره آتی 
درمجموع بر تعداد روزهای همراه با بارش سنگین در اکثر ایستگاه‌های موردمطالعه افزوده خواهد شد. خزانی 
و همکاران (۱۳۹۸) با استفاده از مدل HadGEM2‏ به پیش‌نگری تغییرات دما و بارش در حوضه سیرا 
پرداختند. نتایج نشان داد كه در دوره آینده بارش کاهش و دما افزایش می‌یابد. حیدری و همکاران(۱۳۹۸) 
چشم‌انداز تغییرات بارش را در منطقه شمال غرب با استفاده از مدل‌های گردش عمومی جو و سناریوهای 
RCP‏ مورد ارزیابی قرار دادند. نتایج نشان داده که بیشتر مدل‌ها برای دوره‌های آینده افزایش بارش را پیش - 
بینی می‌کنند. کثیری و همکاران (VE)‏ به پیش‌نگری تغییرات بارش در سواحل جنوبی دریای خزر با 
استفاده از سناریوهای مختلف واداشت تابشی پرداختند. نتایج نشان داد بارش در سواحل جنوبی دریای خزر 
تا پایان قرن ۲۱ نسبت به دورة پایه به‌طور متوسط ۵/۷ درصد افزايش خواهد یافت. 

جمع‌بندی پیشینه تحقیق نشان داد که بیشتر مطالعات انجام شده در زمینه تغییر اقلیم در ایران و در منطقه 
شمال غرب با استفاده از برونداد مدل‌های گزارش سوم (CMIP3)‏ و يا گزارش پنجم (CMIPS)‏ و داده‌های 
اين مدل‌ها صورت گرفته است و بررسی بارش به‌صورت جامع و با استفاده از مدل‌های گزارش ششم 
6 صورت نگرفته است. لذا در اين پژوهش درستی سنجی مدل‌های متداول و پرکاربرد گزارش ششم 
در دسترس با استفاده از متداول‌ترین سنجه‌های آماری بررسی شد و تغییرات بارش در شمال غرب ایران در 
دوره آینده (۲۰۲۱-۲۱۰۰) با استفاده از خروجی مدل منتخب با کمترین میزان خطا تحت سناریوهای مختلف 
SSP‏ بيش نكرى گردید. 

با توجه به آن که منابع آب در معرض خطرات ناشی از تغييرات اقلیم قرار دارند. بررسی تغییرات بارش در 
سال‌های آینده می‌تواند راهگشای معضلاتی همچون خشکسالی. سیلاب‌های ناگهانی و تخریب‌های 
زیست محیطی باشد. 


Y‏ - داده‌ها و روش‌شناسی 
۱-۲- منطقه موردمطالعه 


محدوده موردمطالعه اين پژوهش. شمال غرب ol ul‏ است که از عرض جغرافیایی ۳۶ درجه و ٤٤‏ دقیقه 
الى YA‏ درجه و ۶۷ دقيقه عرض شمالی و طول جغرافیایی ٤٤‏ درجه و ۵ دقیقه الى £A‏ درجه و ۵۲ دقیقه 


داده‌های بارش Yo‏ ايستكاه همدیدی به‌صورت روزانه» با توزيع مناسب و دارای دوره آماری بلندمدت 


سال سیزدهم پیش نگری بارش با استفاده از مدل‌های 0811276 تا SLL‏ قرن ۲۱ در شمال غرب ایران ۱۷ 


(۱۶ ۱۹۸۵-۵۰ حداقل سی‌ساله در سطح منطقه موردمطالعه استفاده LŚ‏ که موقعیت جغرافیایی Lal‏ در 
شکل (۱) آورده شده است. 
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شکل ۱- موقعیت جغرافیایی و تویوگرافی منطقه موردمطالعه و ایستگاه‌های هواشناسی موردبررسی 


-Y-Y‏ مدل جهانی (GCM)‏ و سناریوهای موردبررسی 

در اين پژوهش از داده‌های بارش ایستگاه‌های همدید شمال غرب کشور به‌صورت روزانه در دوره آماری 
تلد و ساسا با کرو قاری ستل‌های جما 548-9313 اقات كلد که ان ماران هر اکتا 
کشور دریافت گردید. جهت بررسی تغییرات بارش در دوره‌های آینده موردبررسی (۲۰۲۱-۲۱۰۰) ابتدا ۷ 
مدل GCM‏ از CMIP6‏ که دارای داده‌های بارش با وضوحبالا بودند اتتخاب شدند. مشخصات مدل‌های 
موردبروسی در جدول (۱) آوزده شده است, داده‌های بارش مدل‌های مذ‌کور از مرکر ۳50 دانلود گردید, 
سپس polie‏ بارش مشاهداتی و تاریخی برای هر یک از ایستگاه‌های موردمطالعه با تهیه پرنامه در محيط 
0 با استفاده از نزدیک‌ترین كريد داده‌های GCMs‏ در دوره پایه (۱۹۸۵-۲۰۱۶) استخراج 
گردید. سپس از روش تصحيح اریبی يا روش مقیاس کاهی خطى جهت ریزگردانی داده‌ها استفاده گردید. در 
ا ی اهداق و RUN E‏ ا ای رر اکن ریات غا 


1 https://esgf-node.lInl.gov/search/cmip6/ 


\VA‏ جغرافیا و مخاطرات محیطی شمارة اول 


(RMSE)‏ ميانكين مطلق خطا (MAE)‏ و ضریب تعیین Rsqrt)‏ ارزیابی گردید. يس از ارزیابی مدل‌های 
مختلف GCM‏ با استفاده از بهترین مدل» تغییرات بارش در دوره آینده (۲۰۲۱-۲۰۵۰) تحت سه سناریوی 
واداشت اجتماعی- اقتصادی (SSP)‏ پیش‌نگری گردید و تغییرات آن در دوره آینده نسبت به دوره پایه 
محاسبه گردید. دلیل انتخاب مدل‌های مذکور از مجموعه مدل‌های CMIP6‏ اعتبار دقت و در دسترس بودن 
سناریوهای موردنیاز جهت مطالعه است. ينج سناریوی SSP‏ از 5581 تا 85۳5 وجود دارد که به ترتیب 
نشان‌دهنده گذر از جاده سبز» وسط جاده جاده سنگی. جاده تقسیم شده و گرفتن بزرگراه است. کمیت‌های 
SSP‏ اساس تلاش مشترک بين IAV‏ و جامعه IAM‏ ساخته شده است. سناریوهای موردمطالعه در اين 
پژوهش SSP126 ge‏ ۹5۳265 و 55۳585 نشان‌دهنده انتهای يايين محدوده مسیرهای اجباری آینده 
انتهای متوسط تا بالا محدوده مسیرهای اجباری otal‏ و انتهای بالای محدوده آینده هستند (زهرائی و 
حسینی, ۱۳۹۹). در ادامه با استفاده از آزمون آماری من کندال و تخمینگر شیب سن, روند و شيب تغییرات 


Pn‏ مدل موسسه/ کشور قدرت تفکیک 

MRI-ESM2-0 ۱‏ موسسه تحقیقات جوی, ژاپن APT‏ 
MPI-ESM1-2-LR Y‏ موسسه ماكس پلانک. آلمان ۰/4/4 
GFDL-ESM4 ۳‏ آزمایشگاه ژئوفیزیک و دینامیک سیالات ایالات‌متحده Yox)‏ 
oe el MIROC6 t‏ ۱ ۱۰۶ 
MIROC-ES2L 0‏ واس رس M‏ و ۷۲/۸۸ 
FGOALS-g3 1‏ موسسه فیزیک om‏ جين ۱۱ 

CanESM5 ۷‏ مرکز تحلیل و مدل‌سازی اقلیمی. کانادا ۱/۱۱ 


منبع: ژانگ و همکاران (۲۰۲۱) 
-tY‏ مقیاس AS‏ 


جهت مقیاس کاهی داده GCMs‏ از روش مقیاس GAS‏ خطى استفاده گردید. مقیاس کاهی بر اساس 
ميانكين تفاوت بين سری‌های زمانى مشاهده شده ماهانه و سری‌های زمانی دوره تاریخی GCM/RCM‏ در 
همان دوره از سرى زمانى مشاهده شده است. سپس اين تفاوت‌ها براى داده‌های آب‌وهوایی شبیه‌سازی شده 


برای دریافت متغیرهای آب‌وهوایی با اریبی اصلاح شده اعمال می شوند (شرستا و همکاران۲۰۱۵). رويكرد 


1 root-mean-square error 
2 mean absolute error 


سال سيزدهم پیش نكرى بارش با استفاده از مدلهاى CMIPO‏ تا يايان قرن 7١‏ در شمال غرب ايران WA‏ 


مقياس بندی خطی (LSBC)‏ با مقادير اصلاح شده ماهانه بر اساس تفاوت بين مقادير مشاهداتى و 
شبیه‌سازی شده عمل می‌کند. داده‌های بارش با ضریبی بر اساس نسبت میانگین بلندمدت ماهانه داده‌های 
مشاهداتى و داده‌های کنترل اصلاح می‌شود. در این حالت شرط بر اين است که ضرایب تصحیح اعمال شده 


وافزودهها حتى برای شرايط آتى نيز بدون تغيير باقی می‌مانند (توجباين و همکاران ۲۰۱۲). 


Um Pops (t) \ 
Pee USP ee cy, | مستت‎ | Q) 


Hm Pcontr (t) 


PRE) = Po [East 00 


Hm Pcontr (t) 


در روابط P YL‏ بارش است؛ contr‏ تعداد سری‌های شبیه‌سازی شده GCM‏ در طول »55« کنترل Obs‏ 
دوره زمانی مشاهداتی fre‏ سرس زمانی پیش‌نگری شده آینده است که بايد اریبی ol‏ حذف شود؛ BC‏ سری 
زمانی پیش‌نگری شده آینده است که اریبی آن حذف شده است. t‏ گام زمانی و Hm‏ ميانكين بلندمدت ماهانه 
است (منذر و (Y Yr col s‏ با اعمال روش مقیاس کاهی بر روی داده‌های شبیه‌سازی شده عملکرد مدل 


در شبیه‌سازی بلندمدت متغیر بارش به مقدار زیادی افزایش می‌پابد. 
۵-۲- آزمون ناپارامتریک آماری من کندال 


آزمون ناپارامتری من- کندال که توسط من در سال ۱۹2۵ ارائه و سپس توسط کندال در سال ۱۹۷۵ 
تکمیل گردید. بر Gal‏ مرتبه داده‌ها در یک سری زمانی استوار است. اين آزمون برای بررسی تصادفی بودن 
داده‌ها pte)‏ وجود روند) در مقابل وجود روند در سری‌های زمانی هیدرولوژیکی و هواشناسی مورد استفاده 
قرار می‌گیرد (ژانگ و همکاران. ۲۰۰۰). Cope‏ اين آزمون نسبت به ple‏ آزمون‌های تعیین روند استفاده از 
مرتبه داده‌ها در سری زمانی بدون در نظر داشتن مقدار متغیرها می‌باشد که به دلیل وجود چنین خاصیتی. 
می‌توان از اين آزمون برای داده‌های دارای چولگی نیز استفاده کرد و نیازی نیست که داده‌ها در قالب توزیع 
حاصی در آیند. فرض موردبررسی در اين آزمون به‌صورت زیر می‌باشد: 

Ho‏ داده‌ها به‌صورت تصادفی توزیع شده‌اند (روند ندارند). 

Hy‏ داده‌ها دارای روند هستند. 


در اين i gal‏ اگر X1,X2,.--5Xn‏ مشاهدات مورد نظر باشند» آنگاه داریم: 


1 Linear Scaling Bias Correction 


۱۸۰ جغرافیا و مخاطرات محیطی شمارة اول 


71-1 n 


uS Y'sm(s, - x,) (v) 


k=l j=k+1 
تابع علامت می‌باشد که به‌صورت زیر تعریف می‌گرده:‎ sgn كه رابطه فوق 1" مقدار داده زام و‎ 
+1 if (x, -x,)»0 
sgn(x)- 40 if (x, -x,)-0 (£) 
-1 dr (x, -x,)«0 


J 


در واقع در اين آزمون هر داده با تمامی داده‌های پس از خود مقایسه می‌شود؛ در این al>‏ می‌توان 
بجای استفاده از مقادیر اصلی دادههاء از مرتبه داده‌ها درمجموعه مورد نظر (سری زمانی) استفاده کرده و 
مرتبه‌ها را به همین روش مقایسه نمود. به اين دلیل. اين آزمون یک آزمون بر پایه مرتبه داده‌ها (Rank base)‏ 
می‌باشد. با فرض این که داده‌ها مستقل بوده و توزيع يكنواخت دارند. ميانكين و واریانس 5 از روابط زیر به 


دست مىآيد: 


E(S) =0 (0) 


n(n -1(2n 4 5)- X r(t, -1)(2t, +5) 5 


Var(S)= T 


در رابطه بالا: 

n‏ تعداد دادههاء m‏ تعداد گره‌ها و t‏ تعداد osla‏ در هر گره می‌باشد. 

منظور از گره اين است که اگر از یک مقدار داده, بیشتر از یکی وجود داشته باشد. این polite‏ مساوی. 
تشکیل یک گره را می‌دهند و تعداد اين polio‏ مساوی در ply elm oS‏ ) می‌باشد. 


آماره اين آزمون (Z)‏ دارای توزیع نرمال بوده و از رابطه زیر به دست می‌آید: 


5-1 
i S>0 
JVar(s) f 
7 - 40 if S=0 (V) 
S+1 s<0 


فرض صفر پذیرفته می‌شود و در غير اين صورت. فرض صفر رد خواهد شد. در حالت رد فرض صفر 


سال سيزدهم پیش نكرى بارش با استفاده از مدلهاى 081156 تا يايان قرن 7١‏ در dius‏ غرب ايران ۱۸۱ 


(وجود روند)» درصورتی که 0 باشد. سری زمانی دارای روند مثبت (صعودی) و درصورتی که 0 باشد. 
|Z| > 1.64‏ در سطح معنی‌داری E‏ درصد. اگر 


zo Z| 2.57 <€‏ . 
در سطح معنی‌داری ۵ درصد و اگر pe‏ | در سطح معنی‌داری ۱ درصد فرض صفر (عدم 


وجود روند) رد می‌شود A(\AV0 JS)‏ 


Ge‏ زمانی دارای روند منفی (نزولی) خواهد بود. اگر 
|Z| > 1.96‏ 


ددسو كيدا 

تخمين شيب خط روند dau‏ سن (OAM)‏ برای تشخيص روند در یک سرى زمانى داده‌ها ارائه شد. اين 
روش یک روش نايارامترى است که با استفاده از تحليل تفاوت بين مشاهدهاى یک سرى زمانی» یک شيب 
ميانه را برای سری زمانی داده‌ها از رابطه زیر محاسبه می‌کند. در روش سن برای هر دو جفت داده متوالی 
یک شیب محاسبه می‌شود. شیب‌های محاسبه شده تشکیل یک سری زمانی می‌دهند که میانه laol‏ شیب خط 
روند را OLY‏ می‌دهد. در آزمون من- کندال شیب خط روند سری داده‌ها با روش تخمین گر سن از رابطه 


ناپارامتری زیر محاسبه می‌گردد. 


"UN ar c PP 

B = Median | ۳ > i) 
J-1 (A) 
ام می‌باشند و ريك‎ j مشاهداتی 1ام‎ polis برآوردگر شیب خط روند 6 زک به ترتیب‎ B كه در آن‎ 


واحد زمانی بعد از i‏ می‌باشد مقادیر مثبت و منفی 8 نشان دهنده روند افزایشی و کاهشی در سری داده- 
هاست (سن» 8 


۷-۲- ارزیابی و صحت سنحی عملکرد مدل‌ها 

يس از استخراج خروجی مدل‌ها برای صحت سنجی آن‌ها از سنجه‌های متداول آماری استفاده شد که در 
ادامه شرح مختصری از اين شاخصها آورده شده است. در رابطه (A)‏ که تحت عنوان ضریب تعیین (Rsqtt)‏ 
شناخته می‌شود معیاری بدون بعد است و بهترین مقدار آن برابر یک می‌باشد (صداقت کردار و همکاران 
(WV ۷‏ مجذور میانگین مربعات خطا (RMSE)‏ به‌عنوان قیاسی برای OUS‏ دادن اختلاف بين مقادیر 
شبیه‌سازی شده از مقادیر اندازه‌گیری به کار می‌رود اين shee‏ که به‌صورت رابطه (۱۰) تعریف می‌شود 


به‌عنوان مرسوم‌ترین شاخص خطا به کار می‌رود (لين و همکاران» (Y V‏ میانگین مطلق خطا (MAE)‏ برای 


1 Sen's Slope Estimator 


\AY‏ جغرافیا و مخاطرات محیطی شمارة اول 


مقایسه عبارت به عبارت خطای نسبی مقادیر شبیه‌سازی شده با توجه به polis‏ اندازه‌گیری شده به کار مى- 


رود که به‌صورت رابطه (۱۱) ارائه م ىكردد (هو و همكاران» Que‏ 


$xx, 
R? 55 i=l 
N 
20, 2K, 
5 (4) 
FN 
مع وح‎ - X.Y 
RMSE = (| -= 
N (V+) 
N 
2 
MAE = id 
N QD 


در روابط فوق Xo‏ داده‌های مشاهداتى» X,‏ داده‌های شبيهسازى Xo cold‏ ميانكين داده‌های مشاهداتی. N‏ 


تعداد داده‌ها می‌باشد. 


۳- نتایج و بحث 


ay gut‏ دست آمده از اعمال آزمون من-کندال و تخعمینگر شیب سن بر روی داده‌های بارش سالائه در 
دوره آماری GL‏ (۱۹۸۵-۲۰۱۶) در ایستگاه‌های موردمطالعه نشان می‌دهد که بارش در ایستگاه‌های اردبیل» 
يارس آباد. خوی» جلفاء زنجان» سراب. سردشت و ماكو دارای روند افزایشی و در سایر ایستگاه‌ها دارای روند 
کاهشی است. بر اساس آزمون من کندال روند افزایشی بارش در هيج كدام از ایستگاه‌های موردمطالعه معنی- 
دار نشده است و تغییرات مذکور از نوع نوسانات کوتاه‌مدت آب‌وهوایی در منطقه می‌باشد؛ اما روند کاهشی 
در ایستگاه‌های سنندج و مراغه در سطح ۰/۰۱ درصد و در ایستگاه‌های سهند و مریوان در سطح ۰/۰۵ درصد 
معنی دار شده است و حاکی از تغييرات کاهشی بارش در ایستگاه‌های مذكور است. بیشترین شیب تغییرات 
کاهشی نيز مربوط به ایستگاه مریوان با شیب سن برابر ۱۹/۲- است. بیشترین شیب تغییرات افزایشی نيز 
مربوط به ایستگاه سردشت با 1/۵ است. درمجموع تغییرات کاهشی بارش در منطقه موردمطالعه بیشتر از 


تغییرات افزایشی است. (جدول AY‏ 


سال سيزدهم بيش نكرى بارش با استفاده از مدلهاى 081186 تا بایان قرن 7١‏ در شمال غرب ايران Wr‏ 


جدول ۲- نتايج حاصل از آزمون‌های آمارى بر روى سرىهاى بارش سالانه ایستگاه‌های موردمطالعه 


رديف ايستكاه ضريب من كندال شيب سن رديف ايستكاه ضريب من كندال شيب سن 
۱ اردبیل M YS WAR!‏ زرینه ۱/۹۲ ۳۱/- 
۲ ارومیه ۰/۲۹۰ -£6/* IE‏ زنجان ۰/۵۶ ۰/۸۵ 
۳ اهر M S VM-‏ سراب YAV Vit‏ 
ob 3‏ -۱/۹۷ -۱/۵ 1۷ سردشت UNY‏ 1/۵ 
۵ يارس آباد Vey JAY‏ ۸ سقز ۱/۸ 4/47- 
1 پیرانشهر ۰/۲۸۰ ۰/۷۲ ۱۹ سنندج -—VYA HRY PY‏ 
۷ تبریز AD‏ ۱/۲۸۰ ۳۰ سهند -۲/۰۳۲ # Y/* Y4-‏ 
۸ تکاب ۱/۳۷ ۲/۸۹ ۲۱ قروه -۰/۷۵ -MYY‏ 
4 جلفا 1/0۸ ۳/۳۰ ۳۲ ماكو A‏ ۳/۹ 
E‏ خرمدره ۰/۲۱۰ UYM‏ ۲۳ مراغه -0/V4 RY /AQ-‏ 
Mw‏ خلخال -۱/۸۵ -۱/۵۲ yé‏ مريوان *YY-‏ ۷- 
Ww‏ خوى ZA iY‏ 
Yo‏ مهاباد -۰/۹۳ -Y‏ 
۱۳ ميانه -۰/1۵ ZA‏ 
#معنی‌داری در سطح ۱ درصد #معنی‌داری در سطح 6 درصد 


توزیع فضایی تغییرات بارش نيز بر اساس آزمون من کندال و تخمینگر شیب سن نشان می‌دهد که به‌جز 
oly‏ شمالی منطقه مورد و منطقه ارسباران و همچنین مناطقی به‌صورت پراکنده در شرق و غرب منطقه در 
yhe‏ مناطق بارش دارای روند کاهشی است که بیشترین تغییرات أن نيز مربوط به نواحی جنوبی و جنوب 
غربی منطقه موردمطالعه و همچنین جنوب شرق دریاچه ارومیه است. 

در اين پژوهش, ۷ مدل GCMs‏ مورد ارزیابی قرار گرفت و پس از اصلاح اریبی به روش 1.5130 به 
cli!‏ عملگرد مد ها يد بروسی etal‏ بارش Sala,‏ و قاری مدلها در در ا (Y8483Y-15)‏ 
پرداخته شد. نمونه‌ای از اصلاح اريبى بر روى دادههاى بارش ايستكاه تبريز در شکل (۳) آورده شده است. 
همانطور كه مشاهده مىشود دادههاى بارش داراى عدم قطعيت زيادى است لذا انجام عمل اصلاح اريبى 
توانسته است به ميزان زيادى عملكرد مدل را در شبيهسازى بلندمدت بارش در منطقه موردمطالعه افزايش 
دهد. با توجه به تعداد زياد ایستگاه‌ها و نمودارهاى حاصل» فقط خروجى ايستكاه تبريز به‌صورت نمودار در 
اين بخش ارائه كرديد. برای مطالعه بيشتر مىتوان به مطالعات شرستا' (۲۰۱۵) و تتجبين' (Y)‏ در اين 


زمينه مراجعه كرد. 


1 Shrestha 
2 Teutschbein 


VAL‏ جغرافيا و مخاطرات محيطى شمارة اول 


— Corrected 


شکا polis -Y‏ مشاهداتی. تاریخی و اصلاح شده بارش در ایستگاه تبریز در دوره ak‏ (۱۹۸۵-۲۰۱۶) 


يس از مقياس كاهى به نتايج حاصل از ارزيابى مدلهاى مختلف با استفاده از شاخص مختلف در 
ایستگاه‌های موردمطالعه پرداخته شد. نتايج نشان داد دقت اكثر مدل‌های موردبررسى در مناطق كم بارش و 
پربارش منطقه موردمطالعه كمتر از مناطق So‏ است و می‌توان بیان كرد تغييرات زياد بارش توي وكرافى و 
سيستمهاى جوى مختلف از مهم ترين دلايل افزايش حطا در oer‏ ایستگاه‌ها و مناطق موردمطالعه است. با 
توانند الكوى بارش را به‌طور دقيق شبیه‌سازی كنند. درمجموع و بر اساس R? polis‏ مدل MPI-ESM1-2-LR‏ 
با داده‌های بارش مشاهداتی در منطقه موردمطالعه ply‏ با ۰/۸۹ و ميزان RMSE‏ أن نيز برابر با ۱۹/۷ است. 
نتایج رتبه‌بندی عملکرد مدل‌های موردبررسی در شبیه‌سازی بارش شمال غرب کشور در جدول Y‏ آورده شده 
است. يس از مدل MPI-ESMI-2-LR‏ دو مدل MIROC-ES2L‏ و MIROC6‏ هر دو با ضریب تعیین برابر با 
۵ بهترین مدل‌های CMIP6‏ برای پیش‌نگری بارش در منطقه موردمطالعه هستند. همچنین مدل‌های 
MIROC-ES2L‏ و FGOALS-g3‏ در حد متوسط هستند. مدل CanESMS‏ نیز با R?‏ برابر با ۰/۶۹ ضعیف- 
ترين عملکرد را در بين مدلهاى موردبررسى جهت شبيهسازى بارش در منطقه موردمطالعه داشته AL‏ بر 
اساس شاخحص RMSE‏ و MAE‏ نيز مدل MPEESMI-2-LR‏ بهترين عملكرد را در بين مدلهاى 
موردبررسی داشته TE aul‏ در oer‏ مناطق دارای کم برآوردی و در برحی مناطق pe‏ دارای بیش 
برآوردی بوده است. بر اساس دو شاخص مذکور نیز مدل FGOALS-g3‏ با RMSE‏ برابر با ۶۰/۳ و MAE‏ 
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جدول ۳- رتبه‌بندی مدل‌های 021176 جهت شبيهسازى بارش بر اساس ميانكين سنجه‌های خطاسنجی 


ردیف مدل ضريب تعيين | مجذور میانگین مربعات خطا | میانگین مطلق خطا 
MPI-ESM1-2-LR ۱‏ كلل ۱۹/۷۱ ۱9/۰ 
MAE YY/VO */A0 MIROC-ES2L Y‏ 
YY/AA ۳۸/۰۵ ۰/۸۵ MIROC6 ۳‏ 
VINE AY MRI-ESM2-0 t‏ ۱۹/۸۹۶ 
YYAA ۰/۲ (MO FGOALS-g3 ۵‏ 
YAAA YVOA "oY GFDL-ESM4 1‏ 
CanESMS M‏ ۹ 10/4۱ ۲۰/۳۹ 


يس از انتخاب بهترین مدل با بهترین عملکرد در شبیه‌سازی بارش منطقه موردمطالعه یعنی مدل MPI-‏ 
ESMI-2-LR‏ به پیش‌نگری بارش با استفاده از مدل مذکور در دوره‌های آینده (۰۲۰۲۱-۲۰۵۰ ۲۰۵۱-۲۰۸۰ 
و ۲۰۸۱-۲۱۰۰) براساس سه سناریوی موردبررسی پرداخته شد. نتایج نشان داد که بر اساس سناریوی 
5581-6 بارش در منطقه موردمطالعه در دوره (۲۰۲۱-۲۰۵۰) بين ۲۰۰ تا ۸۵۰ میلی‌متر خواهد بود که 
نسبت به دوره ah‏ به‌طور متوسط ۳/۲ درصد YN)‏ میلی‌متر) در منطقه افزايش OUS‏ می‌دهد. در دو دوره 
دیگر نیز (۲۰۵۱-۲۰۸۰ و ۲۰۸۱-۲۱۰۰) میزان بارش به‌طور متوسط به ترتیب ۰/۰۳ و ٤/٦‏ درصد نسبت به 
دوره پایه افزايش خواهد یافت که نسبت به دوره قبلی کاهش زیادی را OUS‏ می‌دهد. بر اساس توزیع فضایی 
بارش در دوره‌های آینده در هر سه دوره موردبررسی بیشترین تغییرات افزایشی مربوط به بخش‌هایی شمالی 
منطقه موردمطالعه یعنی منطقه جنگلی‌های ارسباران است که می‌تواند به دلیل پوشش جنگلی منطقه» نزدیکی 
به منابع آبی و همچنین توپوگرافی منطقه باشد. در ساير مناطق موردمطالعه میزان بارش کاهش پیدا خواهد 
کرد که بیشترین میزان آن مربوط به نواحی جنوب غربی منطقه موردمطالعه می‌باشد که جزو مناطق پربارش 
در منطقه موردمطالعه نیز هستند. درمجموع میزان بارش بر اساس سناریوی SSP1-2.6‏ در هر سه دوره آینده 
نسبت به دوره پایه در منطقه موردمطالعه افزایش بيدا خواهد کرد که میزان OF‏ به‌طور متوسط در دوره ۲۰۲۱ 


تا ۲۱۰۰ برابر با ۲/٩‏ درصد خواهد بود (شکل ۵). 
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شکل ۵- توزیع فضایی بارش در منطقه شمال غرب کشور بر اساس سناریوی SSP1-2.6‏ در سه دوره آینده 


نسبت به دوره پایه 


توزیع فضایی مقادیر بارش بر اساس سناریوی 5523-0 نشان داد که بارش در منطقه موردمطالعه در سه 
دوره آینده موردبررسی (۵۲۰۲۱-۲۰۵۰ ۲۰۵۱-۲۰۸۰ و ۲۰۸۱-۲۱۰۰) به ترتیب» 0۲۹/۸ ENA‏ و ۸۲/۳ میلی‌متر 
کاهش خواهد يافت كه نسبت به دوره يايه به ترتيب الى 1/۲ و YY/Y‏ درصد کاهش را نشان می‌دهد. 
کمترین تغییرات کاهشی حاصل مربوط به نواحی شمالی منطقه و بیشترین تغییرات مربوط به نواحی جنوب 
منطقه موردمطالعه کاهش بيدا خواهد کرد که اين میزان کاهش در سطح منطقه در دوره ۲۰۲۱ تا ۲۱۰۰ نسبت 


به دوره al‏ به‌طور متوسط برابر با ۵ درصد خواهد بود (شكل na‏ 
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شکل - توزیع فضایی بارش در منطقه شمال غرب کشور بر اساس سناریوی SSP3-7.0‏ در سه دوره آینده 


نسبت به دوره پایه 


دوره (۲۰۳۱-۲۰۵۰) بين ۰ نا ۰ ميلى متر خواهد بود که نسبت به دوره ak‏ ۸ درصد (۱/۱میلی‌متر) 
افزایش خواهد داشت؛ اما در دو دوره دیگر نيز (۲۱۵۱-۲۰۸۰ و ۰--۲۰۸۱) میزان بارش به‌طور متوسط در 
منطقه موردمطالعه به ترتيب» ۲ و ۶/۵ درصد نسبت به دوره ab‏ کاهش خواهد یافت. بیشترین تغییرات 
افزایشی و کاهشی در پهنه موردمطالعه نیز همچون دو سناریوی AS‏ به ترتیب مربوط به نواحی شمالی و 
c £g‏ غربى منطقه خواهد بود. درمجموع ميزان بارش در يهنه موردمطالعه بر اساس اين سناریو در دوره 
۲۱ تا ۰ نسبت به دوره پایه Y/A‏ درصد کاهش بيدا خواهد کرد (شکل (V.‏ 


\M‏ جغرافیا و مخاطرات محیطی شمارة اول 
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شکل ۷- توزیع فضایی بارش در منطقه شمال غرب کشور بر اساس سناریوی SSP5-8.5‏ در سه دوره آینده 


نسبت به دوره پایه 


درمجموع بر اساس نتایج حاصل از بررسی سناریوهای مختلف به‌جز سناریوی SSPI-2.6‏ روند بارش در 
دوره‌های موردبررسی در منطقه موردمطالعه کاهشی است و بیشتر تغییرات به مناطق پربارش منطقه 
موردمطالعه اختصاص دارد. تغیبرات افزایشی در برخی مناطق موردمطالعه نيز می‌تواند ناشی از نزدیکی به 
منابع آبى و مسیر رطوبت دریای سياه و مدیترانه و یا کوهستانی بودن و پوشش گیاهی و جنگلی Of‏ منطقه 
باشد. مقایسه سناریوهای مختلف نیز نشان می‌دهد که سناریوی SSP3-7.0‏ تغییرات کاهشی بیشتری نسبت 
به So‏ سناریوها OLE‏ می‌دهد و برعکس سناریوی SSPI-2.6‏ تغییرات کمتری را OUS‏ می‌دهد که ناشی از 
ویژگی‌هایی است که هر یک از سناریوها به نمایش می‌گذارند (شکل AA‏ 
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نكرى بارش با اسفاده از مدلهاى CMIPG‏ تا بایان قرن ۲۱ در dies‏ غرب ایران 


۱۸۹ 


۱۹۰ جغرافیا و مخاطرات محیطی شمارة اول 


*- نتیجه گیری 


بيش نكرى تغییرات بارش برای اطلاع از میزان تغییرات آن در آینده و در نظر گرفتن تمهیدات لازم برای 
تعدیل اثرات سوء ناشی از تغیبرات آب‌وهوایی بر منابع آبی و کشاورزی از اهمیت زیادی برخوردار است. 
اهمیت پیش‌نگری بارش زمانی آشکارتر است که در جایگاه‌های جغرافیایی ریزش ناچیز يا قابل‌توجه و يا 
ناگهانی حجم زیادی از بارش رخ دهد. لذا با توجه به اهمیت موضوع. در اين پژوهش به ارزیابی ۷ مدل 
GCMs‏ از CMIP6‏ در شبیه‌سازی بارش و همچنین پیش‌نگری تغییرات بارش در منطقه شمال غرب کشور 
تحت سه سناریوی واداشت اجتماعی- اقتصادی (SSPs)‏ در سه دوره زمانی آینده (۰۲۰۲۱-۲۰۵۰ ۲۰۸۰- 
۱ و ۲۰۸۱-۲۱۰۰) نسبت به دوره ah‏ (۱۹۸۵-۲۰۱۶) پرداخته شد. نتايج حاصل از ارزیابی عملکرد 
مدل‌های GCM‏ موردبررسی با استفاده از شاخص‌های خطاسنجی مختلف نشان داد که مدل MPI-ESM1-2-‏ 
LR‏ در منطقه موردمطالعه از دقت بیشتری نسبت به ساير مدل‌ها جهت شبیه‌سازی بارش برخوردار است؛ 
بنابراین بر با استفاده از داده‌های مدل MPI-ESMI-2-LR‏ و سناریوهای آن اقدام به پیش‌نگری بارش گردید. 
نتایج حاصل از پیش‌نگری تغییرات بارش نيز نشان داد كه بارش بر اساس سناریوی خوش‌بینانه SSP1-2.6‏ تا 
پایان قرن ۲۱ در سطح منطقه حدود ۲/١‏ درصد نسبت به دوره ab‏ افزایش بيدا خواهد کرد. با اين وجود 
به‌جز ols‏ شمالی منطقه موردمطالعه در سایر پهنه‌های موردمطالعه میزان بارش در دوره آینده نسبت به 
دوره پایه کاهش خواهد یافت که بیشترین میزان کاهش آن نيز مربوط به نواحی پر بارش در جنوب غرب 
منطقه موردمطالعه است. به‌جز ایستگاه‌های مریوان» ارومیه و تبریز و زرینه در ساير ایستگاه‌ها افزایش بارش 
معنی‌دار شده است که بیشترین افزايش در سطح معنی‌داری ۰/۱۱ درصد مربوط به ایستگاه اردبیل می‌باشد؛ 
اما بر اساس سناریوهای SSP3-7.0‏ و SSP5-8.5‏ میزان بارش تا OLL‏ قرن ۲۱ به ترتیب ۱۶/۵ و V/V‏ درصد 
کاهش بيدا خواهد کرد و بیشترین کاهش أن نیز مربوط به نواحی جنوب غربی منطقه موردمطالعه است. بر 
اساس سناریوی SSP3-7.0‏ در بیشتر ایستگاه‌های موردمطالعه در منطقه بارش کاهشی و معنی‌دار شده است. 
بیشترین میزان کاهش بارش مربوط به ایستگاه‌های سردشت. خوی و جلفا می‌باشد که در سطح ۰/۰۱ درصد 
معنی‌دار شده است. نتایج حاصل از سناریوی SSPS-8.5‏ تا پایان قرن ۲۱ نيز نشان داد که بارش در منطقه 
به‌جز ایستگاه میانه در تمامی ایستگاه‌ها دارای روند کاهشی است و میزان اين کاهش نسبت به سناریوی قبل 
اندکی بیشتر است. میزان اين کاهش در مناطق غربی و جنوب غربی منطقه موردمطالعه در ایستگاه‌های 
سردشت. پیرانشهر مهاباده يارس آبادء ماکی جلفا و خوی در سطح ۰/۰۱ درصد معنی‌دار شده است. quU‏ 
حاصل با مطالعات مجدی و همکاران (۲۰۲۲) مبنی بر کاهش بارش در بیشتر مناطق خاورمیانه و شمال 
آفریقا بر اساس متوسط ۲۳ مدل GCM‏ مطالعه حیدری و همکاران(۱۳۹۸) بر مبنای بررسی مدل‌های 


گردش کلی جو و استفاده از سناریوهای RCP‏ مبنی بر کاهش بارش در دوره‌های ( ۲۰۵۱-۷۰ و ۰ 
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۲۰۵۰ و همچنین مطالعات صلاحی و همکاران (۱۳۹۵) مبنى بر کاهش ۲ درصدى بارش در دهه‎ (Ae 
دریاجه ارومیه مطابقت دارد. از نتایج حاصل از اين پژوهش می‌توان در مقابله با معضلاتی‎ soul در حوضه‎ 
همچون خشکسالی. سیلاب‌های ناگهانی و تخریب‌های زیست‌محیطی و همچنین مدیریت منابع آب؛‎ 
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